Astronomy & Astrophysics manuscript no. 1 1653 


© ESO 2009 


March 20, 2009 





0^ 
O 
O 
(N 



O 
(N 

^' 

O 



(N 

in 
m 



A self-consistent empirical model atmosphere, abundance and 
stratification analysis of the benchmark roAp star a Circini^ 

O. Kochukhov', D. Shulyak^, and T. Ryabchikova^'^ 

' Department of Physics and Astronomy, Uppsala University, Box 515,751 20, Uppsala, Sweden 

^ Institute of Astronomy, Vienna University, Turkenschanzstrasse 17, 1180 Vienna, Austria 

^ Institute of Astronomy, Russian Academy of Science, Pyatnitskaya 48, 119017 Moscow, Russia 

Received 13 January 2009 / Accepted 15 March 2009 

ABSTRACT 

Context. Chemically peculiar (CP) stars are unique natural laboratories for investigation of the microscopic diffusion processes of 
chemical elements. The element segregation under the influence of gravity and radiation pressure leads to the appearance of strong 
abundance gradients in the atmospheres of CP stars. Consequently, the atmospheric temperature-pressure structure of these objects 
could deviate significantly from the atmospheres of normal stars with homogeneous abundances. 

Aims. In this study we performed a self-consistent, empirical model atmosphere study of the brightest rapidly oscillating Ap star 
a Cir. We account for chemical stratification in the model atmosphere calculations and assess the importance of non-uniformed verti- 
cal element distribution on the model structure, energy distribution and hydrogen line profiles. 

Methods. For chemical stratification analysis we use the DDAFit minimization tool in combination with a magnetic spectrum syn- 
thesis code. The model atmospheres with inhomogeneous vertical distributions of elements are calculated with the LLmodels stellar 
model atmosphere code. 

Results. Based on an iterative procedure of the chemical abundance analysis of 52 ions of 35 elements, stratification modeling of 
4 elements (Si, Ca, Cr and Fe) and subsequent re-calculations of the atmospheric structure, we derived a new model atmosphere of 
a Cir, which is consistent with the inferred atmospheric chemistry of the star. We find T^g = 7500 K, log g = 4. 1 , and demonstrate that 
chemical stratification has a noticeable impact on the model structure and modifies the formation of the hydrogen Balmer lines. At 
the same time, energy distribution appears to be less sensitive to the presence of large abundance gradients. 

Conclusions. Our spectroscopically determined T^g of a Cir agrees with the fundamental effective temperature of this star. This shows 
that temperatures inferred in detailed spectroscopic analyses of cool magnetic CP stars are not affected by a large systematic bias. 

Key words, stars: abundances - stars: atmospheres - stars: chemically peculiar - stars: individual: a Circini 



^ 1. Introduction a few available theoretical diffusion calculations dBabell [19921: 
^ . iLeBlanc & Moninl2004l:rLeBlanc et al .' 2009). At the same time, 
a> . Atmospheres of certain chemically peculiar (CP) B-F type stars rare-earth elements (REE) show completely opposite picture 
; are dynamically stable against strong convective and turbulent ^ith a large overabundance occuiTing in the superficial atmo- 
> motions. In combination with strong, large-scale magnetic fields spheric layers and normal concentration or underabundance in 
and slow rotation, this provides an ideal con dition for the op - the continuum forming regions (Mashonkina et al. 2005, 2008). 
eration of the microscopic particle diffusion ( |Michaud| | l970|) . i^j-gg abundance anomalies and chemical gradients, 
^ The overall, conspicuously non-solar, chemical signahire of the t^e atmospheric structure of CP stars can deviate significandy 
CP-star atmospheres is determined by an interplay between the fj-om that of the normal B-F stars with similar fundamental pa- 
effects of the gravitational setthng and radiative levitation of rameters and homogeneous atmospheres. For example, the pres- 
chemical species. Due to rapid variation of thermodynamic prop- g^^g core-wing anomaly in the Balmer line profiles of cool Ap 
erties in the outer stellar envelope, radiative diff^usion is also able gj^rs indicates that their tem perature stru cture is most likely dif- 
to produce a vertical element separation within the line-forming fg^gnj f^^^ no^j^al stars ( C owlev et aLllMlt iKochukhov et all 
atmospheric region of certain CP stars. The resuking picture 2002). If unaccounted, these effects may lead to significant sys- 
of chemical stratification depends upon the eff^ective tempera- tematic errors in determination of the stellar atmospheric and 
ture and magnetic field strength of the star (Rvabchikova 2008). fundamental parameters 

However, the Hght and iron-peak elements generally follow a j„ ^^-^^ ^^^^^^ advances in the theoretical treat- 
simple step-like behaviour with a large (relative to solar) element jng„j of chemical stratification in stellar atmospheres 
underabundance m the upper atmosphere and an overabundance ,Hui-Bon-Hoa et al. 200a_fc^lanc & St 
m the lower photospheric layers (see, for example, BagnuloetalJ [Xlecian & StifH 120071: ILeBlanc et"dJ 120091) . existing mod- 
L2001; Ryabchikovaetal. 2002, 2005; Kochukhov et al. 2006D. gj, g^^^^t be applied for the routine analysis of CP-star 
A qualitatively similar chemical stratification is predicted by atmospheres due to their complexity, uncertainties in the 
description of relevant physical processes, incomplete atomic 



5'eMi/o#n>jfre^Mejr.sto:0. Kochukhov, e-mail: olegOfysast. uu.se data, and occasional lack of self-consistency between the 
* Based on observations collected at the European Southern model atmosphere structure and vertical abundance profiles. 
Observatory, Paranal, Chile (ESO program 68.D-0254) Instead, empirical stratification studies, which derive element 



2 



O. Kochukhov et al.: A self-consistent empirical analysis of the benchmark roAp star a Circini 



distributions from the observed line profiles, typically employ 
model atmospheres computed with the assumption of a homo- 
geneous element distribution. This leads to an inconsistency 
between the inferred stratification and the adopted atmospheric 
temperature-pressure structure. 

The aim of the present paper is to explore a new, alternative 
approach to the problem of understanding the interplay between 
non-solar, inhomogeneous chemistry and the atmospheric struc- 
ture of CP stars. Similar to previous empirical stratification stud- 
ies, we use high-quality spectroscopic observations to obtain the 
average chemical composition, vertical distribution of some of 
the most important elements, and place constraints on the stellar 
parameters. Then we make a second step of including individ- 
ual abundances as well as chemical stratificatio n in the model 
atmo sphere calculation with our LLmodels code ( S hulvak et alj 
|2004<) . This two-stage procedure is repeated until convergence is 
achieved, providing us with a self-consistent model atmosphere 
and reliable stellar parameters inferred from best available ob- 
servations. 

We selected the bright (V - 3.2), cool magnetic Ap star 
a Circini (HD 128898, HR 5463, HIP 71908) for our investi- 
gation. This star is the brightest of about 40 currently known 
rapidly oscillating Ap (roAp) stars and is one of the most promis- 
ing objects for detailed asteroseismic in vestigation. It pul sates 
in several non-radial oscillation modes dKurtz et anil994l) and 
was recently observed photometrically by the WIRE satellite 
jBruntt 2007) . The t ime-resolved spectroscopic observations of 
pulsa tions in aCir jRvabchikova et al.l 120071; iKochukhov et alj 
I2OO7 I) allowed to probe the vertical mode cross-section and study 
the dynamics of the outer atmospheric layers enriched in REEs. 
These studies infer a very different pulsation behaviour of the 
lines of iron-peak and rare-earth elements, which yields an inde- 
pendent confirmation of the validity of the chemical stratification 
picture outlined above. The ongoing seismic studies of aCir re- 
quire reliable fundamental parameters and detailed knowledge 
of the stellar atmospheric structure. 

A model atmosphere study of aCir may benefit from a 
rich collection of the photometric, spectrophotometric, and high- 
resolution spectroscopic data collected for this star over several 
decades. In particular, the availability of the observed UV and 
optical spectral energy distribution calibrated in absolute units is 
very valuable for testing theoretical model with chemical strati- 
fication included. The absence of such data leaves significant un- 
certainty in the pa rameter determinatio ns of magnetic CP stars 
(see, for example. [Shulvak et al.ll2008l 1200 9). Furthermore, the 
angular di ameter of a Cir wa s recently measured interferomet- 
rically by iBruntt et all (1200 8', hereafter BNC08). This allowed 
direct determination of the stellar radius for the first time for 
any magnetic Ap star. The effective temperature could be sub- 
sequently inferred from the bolometric flux, almost indepen- 
dently from the model atmosphere techniques (however, still 
using theoretically calculated flux for the IR region due to the 
lack of spectrophotometric observations). Thus, a Cir is a unique 
benchmark cool Ap star and an extremely useful object for the 
verification and calibration of various photometric and spectro- 
scopic techniques of deriving basicstellar properties. In fact, the 
Jeff = 7420 ± 170 K found by iBNC08l is lower than all previ- 
ous photometric and spectroscopic temperature estimates, which 
span the range of 600 K. This new fundamental temperature of 
a Cir differs significantly from Teff = 7900 + 200 K suggested for 
this s tar in the earlier model atmosphere analysis (Kupka et al. 
119961) . raising doubts about the 150-200 K T^fi accuracy often 
claimed in the studies of cool Ap stars. 



The chemic al composition of th e atmosphere of ffCir was 
first studied by 'Kupka et alj ( II99 6'). Based on new fundamen- 
tal values of the atmospheric parameters, BNC08 redetermined 
abundances for a number of chemical elements using models 
with solar metallicity and with individual abundances. However, 
they did not carry out a stratification analysis but assumed 
a chemically-homogeneous atmosphere both for the line for- 
matio n and for model a tmosphere calculation. On the other 
hand, iRvabchikova et all (|2008) presented an analysis of the 
Ca stratifica t ion in a Cir based on the model parameters from 
lKupkaetalJ(ll996h and neglecting the influence of chemical gra- 
dients on the model structure. 

In this paper we determine abundances and reconstruct the 
stratification of most important elements in the atmosphere of 
a Cir. We then explore the effect of this non-solar, inhomoge- 
neous atmospheric chemistry on the derivation of stellar param- 
eters. Since the non-uniform element distribution changes the 
model temperature-pressure structure, we implement an itera- 
tive procedure of the stratification and abundance analysis with 
the subsequent re-calculation of the model atmosphere and re- 
determination of stellar parameters. 

This paper is organized as follows. In the next section we 
give an overview of observational data. Then, in Sect. [3] we de- 
scribe the general methods of chemical abundance and stratifi- 
cation analysis, model atmosphere calculation and atomic line 
data used. In Sect. |4]the results of the stratification modeling are 
presented and the influence of chemical gradients on various ob- 
served stellar characteristics is discussed. The summary of our 
investigation and discussion of our results in the context of other 
recent studies of cool Ap stars are presented in Sect. |5] 

2. Observations 

For abundance and stratification analysis we us ed the spec- 
trum of a Cir obtained with the U VES instrument dPekker et al.l 
I2OOO ) at the ESO VLT. The observations were carried out in 
two dichroic modes, using a 0.5" slit for a spectral resolu- 
tion of 80000. The spectra cover the full wavelength region 
of 3040-10400 A except for a few small gaps. Reduction of 
the UVES spectra w as performed with the automatic pipeline 
dBallester et alJl2000l). The reduction p roces s was described in 
more detail by IKochukhov et alj (|2006|) and IRvabchikova et al.l 
(I2008h . 

The observed ultraviolet spectral energy distribution of 
a Cir was obtained from lUE Newly Extracted Spectra (INES) 
archiv43- We used the average of rebinned, high-dispersion spec- 
tra, recorded with a large aperture in the 1150-3100 A wave- 
length window. For the spectral interval 3200-7500 A we used 
low -resolution energ y di stributions of a Cir from the catalogues 
by BurnashevI d 1 9851) and lAIekseeva et all dl996l) . These data are 
calibrated in absolute flux and show a systematic difference of 5- 
10% below A - 4500 A. The agreement is better than 5% for the 
longer wavelengths. 

The spectrophotometry of a Cir was complemented with 
the broad-band photometric measurem ents in th e UBVRI 
(Johnson etal. 1966) and JHKLM ( Gro ote & Kaufm ann 198j) 
systems. Stellar magnitudes were co nverted to absolute flux 
units with the help of calibrations by iBessell et alj dl998l) and 
lEngels et al.l dl98lh . respectively. 

In addition to using the low-resolution spectral energy distri- 
bution as well as metal and hydrogen line profiles in the UVES 

' |http : // sdc ■ laef f . inta . es/ ines/index2 . html | 
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spectra, we compared predictions of our theoretical models with 
the Stromgren and Geneva photometric parameters of a Cir The 
observed Gen eva colors were obtained from the catalogue by 
iRufenerl (11988 ). For the Stromgren photometry, numerous mea- 
surements listed for a Cir in the SIMBAD databas^ yield av- 
erage parameters b - y = 0.135 + 0.013, mi =0.218 + 0.020, 
ci = 0.783 + 0.006, and H/3 = 2.829 + 0.009. 



3. Methods 

3.1. Abundance and stratification analysis 

The classical abundance analysis in the approximation of 
chemically-homogeneous atmosphere was carried out with 
the help of the equivalent width technique. We used the 
WIDTH9 code feu rucz 1993), modified by V. Tsymbal (see 
iRvabchikova et alT2002l) . to estimate abundances of 35 elements 
for a given model atmosphere. In few cases of blended lines or 
lines affected by the hyperfine structure (hfs) theoretical line pro- 
files calculations were performed with the m agnetic spectrum 
synthesis code SynthMag ( Kochuk hov! |2007|). assuming a ho- 
mogeneous radial field structure with the magnetic field mod- 
ulus of (B> = 2.0kG. This value correspon ds to the upp e r limit 
of t he field strength in a Cir estimated by iKupka et al.l d 19961) 
and BNC08. We also tried (B) = l.OkG and found no significant 
effect on the results of abundance or stratification analysis. 

Since it is not possible to account for the Zeeman splitting 
and polarized radiative transfer in the WIDTH9 code, we used 
a pseudo-mi croturbulence to mimic the magnetic intensification 
effects (see Kupk a et"a]lll99 6). This approach is appropriate for 
the equivalent width analysis of magnetic stars with (B) < 3kG. 
The pseudo-microturbulence parameter of = 1-2 ± 0.2 kms"' 
was estimated by matching profiles of the Fei and Pen lines 
computed with the full treatment of polarized radiative transfer 
and zero microturbulence by a non-magnetic spectrum synthe- 
sis calculation. Since a Cir has an inhomog eneous s urface a bun- 
dance structure (iKochukhov & Rvabchikov a 2001; BNC08I) and 
a non-uniform surface magnetic field distribution, results of the 
abundance analysis and the formal value may slightly change 
with the rotation phase. However, the spectral and magnetic vari- 
ation of this star over its 4.46'' rotation period are mino r com- 
pared to other cool Ap stars (e.g.. IKochukhov et al.l2004l) . which 
excludes a large influence of the horizontal inhomogeneities on 
the results of our analysis. 

Stratification modeling was performed for Si, Ca, Cr, and Fe. 
These elements have sufficient number of lines in the observed 
wavelength region and are able to influence the model structure 
and flux distribution by altering the line blanketing. For these 
four stratified elements we adopted the ste p-function approxi- 
mation of th e vertical abundance distribution (IRvabchikova et afl 
I2005L 12008*). The vertical abundance profile was described by 
four parameters: chemical abundance in the upper atmosphere, 
abundance in deep layers, the position of the abundance jump 
and the width of the transition region where chemical abundance 
changes between the two values. All four parameters were opti- 
mized simultaneously with the IDL-based DDAFit least-squares 
fitting procedure (Kochukhov 2007). This code adjusts the verti- 
cal abundance distribution of a given element to achieve the best 
fit to a set of observed spectral line profiles. In these calculations 
we used the SynthMag code and adopted sin / = 12.5 km s ' . 

We note that the assumption of constant abundances in the 
upper atmosphere and in the lower photospheric layers is related 
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to the lack of spectral lines sensitive to the abundance changes 
in these atmospheric regions. Due to substantial rotation of a Cir 
our choice of lines for the stratification analysis is limited to ab- 
sorption features probing a certain line-forming region, which 
does not extend over the entire atmosphere. In particular, we 
cannot probe layers with log T5000 ^ 0. The cores of strong lines 
could potentially provide us with the information about stratifi- 
cation in the uppermost layers, but such lines are not available 
for the majority of elements. We estimated that robust determi- 
nation of stratification using the lines available in the spectrum of 
a Cir is possible for the log T5000 ~ [-3.0, -1-0.5] region. Outside 
this range of optical depths, an unconstrained solution of the ver- 
tical inversion probl em is non-unique due to ill-conditioned na- 
ture of the problem (IKochukhov et alj|2006l) . Thus, to overcome 
this difficulty, we force chemical abundance to a constant value 
in these layers. 

3.2. Atomic line data 

Atomic data for the abundance calcula tions were extracted 
from the recently u pdated VALD database (iPiskunov et al.l 19951 
iKupka et al.ll 19991 see also http://vald.inasan.ru/~vald/ and ref- 
erences therein). For Tb m and Dy m transition probabilities were 
kindly provided by A. Ryabtsev (private communication). These 
calculations are based on the extended analysis of the energy 
level structure. For several ions (Mni,ii, Coi, Euii, Hon) with 
substantial hfs splitting we included this effect in the abundance 
calculations. The correspondi ng hfs c onstants were taken from 
Kur ucz (2008) f or Mn i, from lHolt et^ (1999.) fo r Mn 11, from 
Pic kerin d(fl996h for Co i from lLawler et all (1200 ll) for Eu 11, and 
from lLawler et al.l (l2004l ) for Ho 11. 

Stratification analysis of Si, Ca, Cr, and Fe requires accurate 
and homogeneous line p arameters for these e lements. For Ca we 
used the same lines as in IRvabchikova et aTl (12008). The line list 
for other elements is given in Table[T] The transition probabilities 
were improved by fitting the solar flux spectrum for a few Si i 
lines and one Fe i line. For Si 11 lines the Stark constants are taken 
from Wilke (2003) , while for al l other lines they are taken from 
Kurucz database (iKuruczll2008l) . 



3.3. Model atmosphere calculations 

Our model atmosphere calculations are based on the ver- 
sion 8.6 of the 1-D, LTE model atmosphere code LLmodels 
dShulvak et al.l |2004|) . This program treats the bound-bound 
opacity by direct, line-by-line spectrum synthesis and is able to 
account for the effects of individual non-solar abundance and in- 
homogeneous vertical distribution of elements. Chemical abun- 
dances and stratification are provided as input parameters for the 
LLmodels code and kept constant in the model atmosphere cal- 
culation process. 

For all model atmospheres calculations in this paper we used 
a 72-layer vertical grid with equidistant spacing in the log T5000 
range between -7 and 2. The frequency-dependent quantities 
were calculated between 50 A and 100000 A with a constant 
wavelength step of 0.1 A. The VALD database was used as 
the main source of the atomic data for compu tation of the Une 
absorption c oefficients. As it was sh own by Kochukhov et al.l 
(I2OO5I) and iKhan & Shulva^ (l2006h . magnetic fields with a 
strength of 1-3 kG expected for aCir could be safely neglected 
in the model atmosphere calculations. Consequently, we ignored 
the effects of Zeeman splitting and polarized radiative trans- 
fer on the atmospheric structure. However, we used a pseudo- 
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Table 1. List of Si, Cr, and Fe spectral lines used for stratification calculations. Columns give the ion identification, central wave- 
length, excitation potential, oscillator strength (log(^/)) and the Stark damping constant (log yst) for T - 10 000 K. The last column 
gives reference for the adopted oscillator strength. 
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microturbulent velocity of 1 kms"' to qualitatively treat mag- 
netic intensification and for consistency with the abundance 
analysis. 

The studies bv lMashonkina et all (l2008ll2005h demonstrated 
that the line formation of rare-earth elements Pr and Nd can 
strongly deviate from the local thermodynamic equilibrium. 
These elements are overabundant in the atmospheres of CP stars 
by several dex compared to the solar atmosphere and have rich 
spectra dominated by the lines of the first and second ions. The 
doubly ionized lines of Pr, Nd, and of some other REEs are un- 
usually strong due to combined effect of stratification of these 
elements and departures from LTE. Therefore, a NLTE stratifi- 
cation treatment of REEs is essential for detailed spectrum syn- 
thesis and also for realistic estimate of their contribution to the 
line opacity in the model atmosphere calculations. However, a 
NLTE analysis is beyond the scope of this paper and cuiTently 
could not be coupled to a model calculation. The large Vg sin / 
of aCir compared to other cool Ap stars also makes it diffi- 
cult to model stratification of REEs since some useful weak 
lines are not accessi ble for measurement. Therefore, we fol- 
lowed IShulvak et al.l (l2009l) in using a simplified treatment of 
the REE absorption. One can approximately reproduce the ob- 



served relative line strengths of the Pr, Nd, Dy, and Tb ions 
by using abundance derived from the second ions and reduc- 
ing the oscillator strengths for the singly ionized REE lines. The 
adopted reduction factor of 1 .5 dex is close to the mean differ- 
ence between concentrations of singly and doubly ionized REEs 
estimated in our classical abundance analysis of a Cir With this 
scaling procedure the Pr, Nd, Dy, and Tb ionization equilibria 
are conserved in the model atmosphere calculations while the 
line strength roughly matches the observed one for both ioniza- 
tion stages. 

3.4. Self-consistent atmospheric modeling 

Previous chemical stratification studies used high-resolution ob- 
servations of CP stars to perform average abundance determina- 
tion and to restore vertical abundance profiles for a few chemi- 
cal elements. Generally, this analysis is performed by interpret- 
ing observed profiles of spectral lines with different strength, 
excitation potential, and ionization degree for a given element. 
However, the shapes and strengths of the line profiles predicted 
by theoretical models may also indirectly depend on the verti- 
cal distribution of chemical elements through the influence of 
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line opacities, modified by peculiar abundances and chemical 
stratification, on the temperature-pressure structure of the model 
atmosphere. Thus, the stratification and abundance analyses are 
linked with the model atmosphere calculation. The extent of the 
feedback of stratification and non-solar chemical composition on 
the model structure cannot be estimated a priori. Consequently, 
one has to apply an iterative procedure of stratification and abun- 
dance analysis followed by re-calculation of the model atmo- 
sphere structure. For the present analysis of a Cir this procedure 
included the following main steps: 

1. calculation of a model atmosphere grid using mean abun- 
dances derived in previous studies; 

2. determination of T^ff and \ogg with the help of spectropho- 
tometry, Balmer line profiles and broad-band colors; 

3. determination of the average abundances and chemical strat- 
ification using high-resolution spectra; 

4. calculation of a new model atmosphere grid, taking into ac- 
count individual abundances and stratification found in the 
previous step; 

5. repeating the entire process starting from step 2 until conver- 
gence is achieved. 

The model atmosphere grids were calculated for the Teff range 
of 7000-8000 K with a 200 K step and for log g = 3.8-4.2 with 
a 0. 1 dex step. We adopted the convergence criterium of the ab- 
sence of significant changes of the mean abundances, stratifica- 
tion profiles of studied chemical elements and model atmosphere 
parameters. 

The final self-consistent, chemically-stratified model atmo- 
sphere is expected to reproduce simultaneously the observed 
photometry, energy distribution, hydrogen line profiles and 
metallic line spectra. The method of self-consistent analysis of a 
cool Ap-star atmosphere applied here is similar to the procedure 
described by Shulyak et al. (2009) for the roAp star HD 24712. 
However, in practice robust re-determination of the atmospheric 
parameters is not possible for the latter star due to the lack of 
observational data. In contrast, a Cir lends itself to a thorough 
study of the interplay between atmospheric chemistry and stel- 
lar parameters. 



: 










L 

: 








- 










J 

















Iteration 










1 








Iteration 


2 






/ 


Iteration 


3 






/ 


Iterotlon 


4 






/ 



















-6 -5 -4 -3-2-10 1 2 




-17 H ....... -i 

-6 -5 -4 -3-2-10 1 2 



-4 " 




-6 -5 -4 -3-2-10 1 2 




-^^ l i 

-6 -5 -4 -3-2-10 1 2 

log T5000 

Fig. 5. Vertical distribution of Si, Ca, Cr, and Fe abundance in the 
atmosphere of a Cir for different iterations of the stratification 
analysis. The thick solid line corresponds to the final element 
distribution. The light horizontal bar shows solar abundances. 



4. Results 

4.1. Chemical abundance pattern 

We started our analysis calculating model atmospheres for a Cir 
with the homogeneously distributed abundances taken from 
[BNC08 and supplemented where necessary by the abundances 
from lKupka et al.l (11996 ). Then the iterative procedure of strati- 
fication analysis was performed following the scheme outlined 
in Sect. 13.41 For each iteration we re-calculated atmospheric 
abundances of all available elements and vertical stratification 
profiles of Si, Ca, Cr, and Fe. In total, four such iterations 
of the abundance analysis and model atmosphere calculations 
were necessary to obtain a converged solution. At each itera- 
tion an adjustment of T^ff by <200 K was required to fit the 
observed spectral energy distribution and the hydrogen line pro- 
files. The final model atmosphere parameters inferred for o-Cir 
are T^^ = 7500 K and log g = 4. 1 . 

A comparison of the observed and calculated profiles of 17 
Fe I and 14 Fe 11 lines used for stratification analysis is presented 
in Fig. [1] These synthetic spectra were computed with the ver- 
tical abundance distribution and atmospheric model obtained at 
the last iteration. Similar plots for Si, Ca, and Cr are provided in 
Figs. |2]^ (Online material). 



The variation of the stratification profiles of chemical ele- 
ments during the iterative process is illustrated in Fig.|5] Neither 
the number of spectral lines used for the stratification analysis 
with DDAFiT nor their atomic parameters were modified during 
the iteration procedure and, thus, all the changes in the stratifi- 
cation profiles from one iteration to another are entirely due to 
changes of the model structure and a small adjustment of Teff. 
In general, the difference between the shape of the stratifica- 
tion profiles obtained on the first and last iterations is relatively 
small for all considered elements. The positions of the abun- 
dance jumps did not change much, although their amplitudes 
increased noticeably for Si and Ca. On the other hand, a large 
change of the stratification parameters occurred at the second it- 
eration, after stratified abundances were introduced in the model 
atmosphere calculation for the first time. For this iteration we 
also had to adopt a 50 K lower effective temperature compared 
to Teff = 7500 K, which provided the best description of hydro- 
gen lines and spectral energy distribution for other iterations. 

For all four elements we obtained a familiar picture of a step- 
like decrease of abundance with height. A slight distortion of 
abundance profiles in deeper layers apparent in Fig. |5]is caused 
by a non-linear mapping from the column mass scale employed 
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Fig. 1. Comparison of the observed Fe line profiles (symbols) and calculations for the final stratified abundance distribution (solid 
line) and for the best-fit homogeneous abundance of Fe (dashed line). 



for vertical inversion to the standard optical depth used in the fig- 
ure. For Cr and Fe the transition region is narrow and is located 
between log T5000 - -0.5 and 0.0. It is centered at log rsooo ~ -1 
and is more extended for Si and Ca. All studied elements show 
a large overabundance with respect to the sun in deep layers 
and significantly (Si, Ca) or marginally (Cr, Fe) sub-solar ele- 
ment concentrations in the upper atmosphere. The formal error 
of the stratification parameters ranges from 0.4-0.8 dex for Ca 
and Si to 0.05-0.15 dex for Fe and Cr, with the upper atmo- 
spheric abundance being more uncertain. At the same time, po- 
sition of abundance jump was determined with the precision of 
better than 0.1 dex in the log T5000 scale for all four elements. 

Table|2] summarizes the final mean abundances derived af- 
ter the fourth iteration of the model atmosphere calculations. 
Element concentrations are reported for 52 ions of 35 elements. 
For completeness we also included mean abundances of the 
stratified element s. Our abundance es timat es are co mpared with 
those reported bv lKupka et alj (1 19961) and lBNCOS' and with the 
solar abundances compiled bv lAsplund e t al. (2005). We find 
that a Cir has a moderate CNO deficiency, almost solar Ca to Fe 
abundances and an overabundance of Cr by « 0.8 dex. The same 
overabundance is obtained for Co when hfs effects are taken into 
account. Ni exhibits a 0.8 dex deficiency. In general , our mean 
abundance results agree well with those by IBNC08I except for 
C, S, Fe, Yii, Zr, and Ba. For the latter the difference in the 
adopted mic roturbule nce parameter might be the reason for the 
discrepancy. |BNC01 determined abundances for three out of 14 
REE elements with stable isotopes, while we succeeded in es- 
timating abundances for 12 REE species. For 6 REEs (Ce, Pr, 
Nd, Tb, Dy, Er) abundances were obtained from the lines of the 
first and second ions separately. Although a Cir has a moder- 
ate REE enhancement of 0.7-1.6 dex according to the results 
obtained using the first ions, the REE anomaly typical for roAp 
stars (Rvabchikova et al. 2004) is clearly observed for all 6 REEs 
where abundances were obtained from the lines of two ionization 
stages. This anomaly is 1.7-1.8 dex for Ce, Pr, Nd, and Tb and 
increases to 2.3 dex for Dy and Er This ionization imbalance 
is an indication of the accumulation of the REEs in the upper 
atmospheric layers (Mashonkina et al. 2005, 2008). Its effect on 
the model atmospher e calculations was investigated in detail by 
IShulvaket al.](l2009h . 



The fourth column in Table|2] gives abundances and stan- 
dard deviations obtained with the WIDTH9 code from n spectral 
lines neglecting hfs effects, while in the next column one finds 
abundances calculated with SynthMag taking hfs effects into ac- 
count. The largest hfs effect is evident for Ho 11 lines, where the 
final abundance decreases by about an order of magnitude. The 
hfs effect on Eu 11 lines is the smallest one and does not exceed 
0.25 dex. The detailed treatment of the hyperfine splitting of the 
Co I lines increases both the discrepancy with Co 11 and the stan- 
dard deviation compared to the WIDTH9 analysis. This could 
be a signature of Co stratification in the atmosphere of a Cir 
Separate hfs analysis of the three moderately strong lines in the 
3500-3600 A region (Wa - 40-82 mA) and five weak to medium 
strength lines in 5300-5400 A region (Wa - 13^2 mA) gives 
log(Co/Mot) = -6.77 + 0.16 and log(Co/Mot) = -6.10 + 0.09, 
respectively. This indicates possible Co accumulation in deep 
layers of the stellar atmosphere. The absence of hfs constants for 
Co II lines does not allow us to carry out a formal Co stratifica- 
tion analysis. In addition, our abundance analysis shows that Sr 
might be strongly stratified in the atmosphere of a C ir similar to 
other Ap stars like, for example, HD 133792 ( Kochukhov et al.l 
I2006h . However, the lack of significant number of suitable un- 
blended lines again prevented detailed stratification analysis for 
this element. 



4.2. Spectral energy distribution and photometry 

We used the observed spectral energy distribution of a Cir as the 
main effective temperature indicator. Previou s attempts to model 
spectrophotometry of magn etic Ap stars dAdelman & Ravi j 
l2000l : lLipski & St^pierill20(^ could utilize only the information 
about the shape of the energy distribution because observations 
were often given in relative flux units. Furthermore, the stellar 
angular diameters were generally unknown for Ap stars and had 
to be treated as free parameters. For a Cir we are in the unique 
position to obtain a very precise estimate of the model atmo- 
sphere parameters from the observed stellar flux thanks to the 
recent direct, interferometric determination of the angular diam- 
eter, e = 1.105 + 0.037 mas (BNC08). Fig.|6]presents a compar- 
ison of the photometric and spectrophotometric measurements 
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Table 2. Abundances in the atmosphere o f aCir with t he err or es timates b ased on n measured lines. Ab undances derived in o ur 
study are compared with those reported bv lKupka et alj (Il996h and lBNCOSl The solar abundances are from lAsplund et afl ( l2005h . 
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of a Cir, discussed in detail in Sect. |2] with the theoretical flux 
distribution computed for the models with individual abundance 
and stratification included. We find that observations are best fit- 
ted by the model with Tetf = 7500 K. The uncertainty of this ef- 
fective temperature estimate, related to fitting ambiguities and 



discrepancy between the two spectrophotometric scans available 
for a Cir, does not exceed 50 K. At the same time, a T^s error 
due to the angular diameter uncertainty is estimated to be 120 K. 
Thus, our analysis of the spectral energy distribution yields an 
effective temperature with the overall accuracy of about 130 K. 
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Fig. 6. Comparison between the observed (symbols) and theoretical (lines) spectral energy distribution of aCir. Photometric and 
spectrophotometric observations, covering the wavelength range from the ultraviolet and optical (left panel) to the infrared (right 
panel), are discussed in the text (see Sect. |2|. The solid line shows theoretical energy distribution computed for the final self- 
consistent model atmosphere of a Cir with parameters Teff = 7500 K and log g = 4.1. The dotted lines illustrate the effect of changing 
Jeff by +200 K. The LLmodels flux spectra are convolved with a FWHM = 10 A Gaussian profile. 
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Fig. 8. Dependence of the synthetic Stromgren 
(left panel) and Geneva (right panel) photo- 
metric parameters on the T^s and log g for the 
models with stratified chemical composition. 
The theoretical temperature {b - y, H/3 for the 
Stromgren and B2 -G, B2- VI for the Geneva 
photometry) and gravity (c[ for the Stromgren 
and d for the Geneva photometry) indicators are 
compared with the corresponding observed pa- 
rameters. Those are indicated by the horizontal 
dashed lines showing ±lcr ranges and by the 
dash-dotted lines (Stromgren photometry only) 
showing the minimum and maximum observed 
values. The vertical dashed lines correspond to 
the 7500 ± 130 K r^ff range. 



Following iBNCOSl we adopted log g = 4.1 for the modeling 
of the observed flux of a Cir. Changes of the surface gravity by 
0.1-0.2 dex lead to a negligible modification of the theoretical 
spectral energy distribution and do not alter results of the T^ff 
determination. 



In the case of a Cir, spectrophotometry is not very sensitive 
to the adopted atmospheric chemical composition and stratifi- 
cation. In Fig. I2] (Online material) we compare LLmodels fluxes 
calculated for our self-consistent, chemically-stratified model at- 
mosphere and for the models adopting homogeneo us solar ab un- 
dance and the chemical composition published by iBNCOSl . All 
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three models predict similar flux distributions, although calcu- 
lations with the non-solar abundances exhibit a more substan- 
tial redistribution of flux from the UV region to the optical and 
IR. The chemically-stratified model predicts a marginally higher 
flux for A < 200 nm compared to the model with homogeneous 
individual abundances. 

Observations in the Stromgren and Geneva photometric sys- 
tems are often used for deter mination of the atmospheric param- 
eters of magnetic CP stars (IHauck & Northlll993l: lNetopil et"aLl 
I2OO8 ). Therefore, we examined an agreement of the photomet- 
ric parameters of a Cir predicted by our self-consistent stratifi- 
cation models with the observed photometric properties of this 
star. For stars with a temperature close to that of a Cir, the b - y 
and HjS parameters of the Stromgren photometry can be used 
as Teff indicators, while ci is mainly sensitive to the su rface 
gravity dMoon & Dworetskvl[T985l: iNapiwotzki et al.lll993l) . For 
the Geneva photometry, the elTective temperature calibrations 
are usually based upon the B2 - G and B2 - VI indices and 
gravity can b e determined with the help of the redden ing-free 
^/-parameter (iHauck & Northlll993t iKunzli et alJll997h . Fig. [8] 
presents these T^g and log g indicators for the entire model at- 
mosphere grid computed with the final abundances and chem- 
ical stratification. The observed range of the Stromgren photo- 
metric observations of a Cir is indicated by the horizontal lines. 
Uncertainties of the Geneva photometric parameters are calcu- 
lated taki ng into account an observational error of 0.004 mag 
(lRufeneii[l988h . 

Generally, Fig. [8] shows an excellent agreement between the 
observed and synthetic b - y, B2 - G, and 82 - VI indices 
for Teff = 7500 + 130 K estimated from the spectrophotometry. 
The Hfi parameter indicates a somewhat lower T^ff on average, 
but it is still consistent with our model predictions if the scat- 
ter of photometric measurements is taken into account. A direct 
spectroscopic comparison between the observed and computed 
hydrogen line profiles presented below provides a more robust 
method to estimate T^g from the Balmer lines. The observed ci 
index suggests log ^ = 4.0 + 0.15 while log ^ = 3.9 + 0.15 fits the 
observed Geneva t/-parameter Both estimates do not contradict 
log g - 4.08 ± 0.09 determined bv BNC08 with stellar evolution- 
ary models. 
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Fig. 9. Comparison of the observed and theoretical Ho' and HjS 
line profiles of a Cir. Observations are shown by the thin line. 
Thick line gives the hydrogen line profiles computed for our fi- 
nal self-consistent, chemically-stratified model atmosphere with 
TtS - 7500 K, log g-W. Dashed lines show the effect of chang- 
ing Teff by +200 K. 



4.3. Hydrogen lines 

Due to its relatively low temperature, the hydrogen line pro- 
files of a Cir are mostly sensitive to T^ff. Thus, we have used 
the Ha and H/3 lines to verify the effective temperature of a Cir 
determined from the stellar flux distribution. The comparison be- 
tween the observed hydrogen lines and the SynthMag calcula- 
tions for the final chemically-stratified models is illustrated in 
Fig-ID An eff'ective temperature of 7470 + 100 K yields the best 
fit to the Ha and Hfi profiles. The error takes into account fitting 
uncertainty and ^1% ambiguity of the continuum normalization 
in the regions around broad hydrogen line wings. 

Unlike the spectral energy distribution, hydrogen lines are 
significantly modified when chemical stratification is included 
in the model atmosphere calculations. Fig. [10] (Online material) 
shows that theoretical Balmer line profiles become systemati- 
cally weaker for the models which neglect chemical stratifica- 
tion. Consequently, the stellar effective temperature will be over- 
estimated by 200 K if such models are to be used for fitting 
the hydrogen line wings. This interesting result is explained in 
Fig.[TT] which shows temperature as a function of continuum op- 
tical depth for models with the same atmospheric parameters but 
different chemical compositions. The calculations with the ho- 



mogeneous individual abundance do not deviate strongly from 
the solar-composition model. On the other hand, the presence 
of a large Cr, Si, and Fe overabundance in deep photospheric 
layers of the chemically-stratified model leads to a significant 
backwarming effect, corresponding to a temperature increase by 
200-400 K for log rsooo ^ 0. The strengthening of the hydrogen 
line wings reflects this change in the model atmosphere struc- 
ture. 



5. Conclusions and discussion 

An empirical investigation of the vertical distribution of chem- 
ical elements in the atmospheres of CP stars is of great impor- 
tance for testing predictions of the microscopic particle diffusion 
theory. Such investigations yield stratification profiles for differ- 
ent chemical elements that could be directly confronted with the 
diffusion calculations. 

In this paper we determined empirical vertical distributions 
of Si, Ca, Cr, and Fe in the atmosphere of the brightest rapidly 
oscillating Ap star a Cir Stratification analysis was comple- 
mented with the most detailed study of average abundances ever 
conducted for a Cir Unlike many previous empirical stratifi- 
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Fig. 11. Temperature as a function of standard optical depth, 
T5000, for the final self-consistent model atmosphere of aCir 
with Teff = 7500 K, log g = 4.1 (solid line) compared to the tem- 
perature structure of the models with the same atmospheric pa- 
rameters but using homogeneous solar (dotted line) and BNC08 
(dashed line) abundances. The upper panel shows the tempera- 
ture differences with respect to the solar-abundance case. 



cation studies of magnetic CP stars, we made an effort to at- 
tain a self-consistency between the inferred chemical proper- 
ties of the stellar atmosphere and its temperature-pressure struc- 
ture. This was achieved by using an iterative procedure of deriv- 
ing stratification and mean abundances from the observed high- 
resolution spectra and subsequent recomputing the model atmo- 
sphere structure. 

Despite a significant improvement of this stellar atmosphere 
modeling methodology compared to previous studies of cool Ap 
stars, the analysis presented in this paper has several limitations, 
stemming from its empirical character The lines available in the 
spectra of aCir can be used to establish stratification of only 
four chemical elements and do not allow exploring stratification 
outside the optical depth range of approximately log tsooo = -3 
to -1-0.5. Furthermore, observations do not allow resolving pos- 
sible variation of chemical stratification profiles over the stellar 
surface, which might be present due to magnetic field effects. 
Potential NLTE effects are also neglected. 

The availability of the photometric measurements in differ- 
ent systems, spectrophotometric scans, hydrogen line profiles 
and, uniquely for a Cir, direct interferometric measurement of 
the angular diameter allowed us to derive accurate atmospheric 
parameters and test the influence of non-solar, inhomogeneous 
abundances on fitting different types of observations. 

The final atmospheric parameters of a Cir were found to 
be Teff = 7500 + 130K and logg = 4.10 + 0.15. Together 
with the chemical stratification and mean abundances derived 
in our study, these parameters reproduce all photometric, spec- 
trophotometric, and spectroscopic constraints available for a Cir 
Our atmospheric parameters are in excellent agreemen t with 
Teff = 7420 + 170 K, log g = 4.09 ± 0.08 determined by BNCOSi 
from the angular diameter measurement, bolometric flux and 



stellar structure calculations. This is an indication that modern 
model atmosphere analyses of cool Ap stars are trustworthy and 
are not affected by significant systematic biases. 

The vertical distribution of chemical elements in the atmo- 
sphere of a Cir was found to be similar to the abundance pro- 
files inferred for other cool Ap stars. All four elements analysed 
for a Cir exhibit a large overabundance in the lower atmospheric 
layers and a solar or subsolar abundance above logTsooo ~ -1. 
We found that the stratification profiles of chemical elements de- 
rived using a model atmosphere computed with homogeneous 
abundances could slightly differ from those obtained using a 
self-consistent analysis. However, in general, these changes do 
not alter significantly the overall energy balance in the stellar 
atmosphere. 

For the atmospheric parameters of a Cir, chemical stratifica- 
tion has a minor influence on the spectral energy distribution in 
comparison to chemically-homogeneous models. Thus, interpre- 
tation of the spectrophotometric scans will not be significantly 
in error if stratification is ignored. Similar result was obtained in 
our study of another cool Ap star HD 24712 (Teff = 7250 K, 
logg = 4.1 Shulvaketal. 2009). We caution, however, that 
the validity of this finding should not be extended outside the 
Tes = 7000-8000 K range addressed in our investigations. Hotter 
Ap stars might have a different abundance stratification that 
makes a larger impact on the atmospheric structure and stellar 
parameter determination. 

On the other hand, we found that the backwarming effect 
introduced by concentration of iron-peak elements in the lower 
atmospheric layers of the chemically-stratified model leads to 
significant strengthening of the hydrogen line wings. The cor- 
responding reduction of Teff inferred from the Balmer lines 
amounts to x 200 K. In the case of a Cir, only chemically- 
stratified models are able to reproduce the Ha and H/3 profiles 
with a Teff consistent with the photometric and spectrophoto- 
metric observations. 

The effective temperature determined in our study is signifi- 
cantly lower than Teff = 790 0+200 K found in the previous spec- 
troscopic analysis of a Cir (Kupk a et al.ll99 6V The origin of this 
discrepancy is twofold. On the one hand, Kupka etal. (199^ 
analysed the star using standard model atmospheres, which ne- 
glected a non-solar chemical composition and element stratifi- 
cation. This led to an overestimation of T^ff obtained with the 
help of photometric techniques and hydrogen line wing fitting. 
On the other hand, the final T^fi was obtained by Kupka et al] 
( Il996h by reducing abundance scatter for a set of Fe i and Fe 11 
lines, assuming a homogeneous distribution of this element with 
height which seemed to be a reasonable approximation at that 
time. We now know that this assumption is incorrect and that, 
due to their different formation heights, the Fe lines of different 
excitation potentials and ionization degree give different appar- 
ent abundance due to stratification effects in cool Ap stars. In 
this situation the criterium of minimum abundance scatter does 
not allow one to recover the true stellar parameters. 

Using the fundamental T^ff and a luminosity of L = 1 0.5 1 + 
O.6OL0 estimated from the Hipparcos parallax, IBNC08I deter- 
mined a stellar mass, M ^1.7 Mq, with the help of stellar struc- 
ture mode ls. We con firmed this result using Padova evolutionary 
tracks dGirardi et alj|200 0). Adopting our spectroscopic T^ff and 
the same luminosity as in .BNC08 . we found M - 1 .69 + O.O2M0 
for the Z = 0.019 grids. The small formal uncertainty of the 
stellar mass is unrealistic given the lack of knowledge about 
the overall metal content of a Cir For example, assuming Z = 
0.008, which is closer t o the recent determin ation of the solar 
atmospheric metallicity jAsplund et al.ll2005l) . yields a mass of 
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M - 1.52 + O.O2M0. For both metallicities the stellar radius 
predicted by the Padova models is in excellent agreement with 
the value of R = 1.967 ± O.O66/?0 determined by BNCOl The 
Z = 0.008-0.019 metallicity interval represents a reasonable es- 
timate o£_th£2angeof_the bulk metal content of nearby field 
stars jBagnulo et al.ll2006h . Unknown initial helium abundance 
and variation in the treatment of convection in stellar structure 
calculations may introduce ad ditional uncertaint ies in the mass 
determination for late-A stars dCunha et alj|2003l) . Thus, despite 
availability of precise estimates of the atmospheric parameters, 
radius and luminosity, we can only conclude that the mass of 
a Cir is likely to be in the range from about 1 .5 Mq to 1 .7 Mq. 

In the light of new, reliable determination of the atmo- 
spheric and fundamental parameters of a Cir, one can use this 
star to test photometric temperature calibrations of cool Ap 
stars. Thes e calibrations were recently summarized and im- 
proved by Netop il et al.l (|2008). For magnetic Ap stars with 
Teff < 9000 K they provided a new B - V calibration and 
recommended using the B2 — G calib ration by iHauck & Norfcl 



1993h 



119821) and the uvby/3 calibration by Na piwotzki et al.l 
for the Geneva and Stromgren photometric systems, respec- 
tively. Application of these photometric temperature calibra- 
tions gives T^siB - y) = 7940 K, Teff (Geneva) = 7678 K, and 
JefrCStromgren) - 7692 K. I n comparison, an ear l ier Str omgren 
photometry calibration by iMoon & Dworetskvl ( Il985h yields 
Teff - 7955 K. Thus, all calibrations overestimate effective tem- 
perature of a Cir by 20 0-450 K while per forming much better 
for other cool Ap stars (iNetopil et al.ll2008D . An unusually high 
abundance of Fe in the atmosphere of a Ci r compared to other 
cool magnetic Ap stars with Teg < 7500 K ('Rvabchikov all2005h 
might be the reason for deviating photometric behaviour of this 
star. 

A large difference in the mean iron-peak element abundances 
of a Cir and of HD 24712, which we analysed using a simi- 
lar s elf-consistent model atmosphere technique (IShulvak et alj 
12009 ). explains the different influence of stratified abundances 
on the model structure of the two stars. Both Fe and Cr are sig- 
nificantly more abundant in a Cir. Consequently, their stratifica- 
tion changes the atmospheric structure to a much greater extent 
than for HD 24712 inspite of a qualitatively similar stratification 
profiles inferred for the two stars. 

Some of the properties of a Cir might be connected with the 
special aspect angle from which this star is visible from Earth. 
The presence of narrow line profiles for the 4.46'' rotation period 
of ffCir indicates a relatively low inclination angle (/ 35°). 
The restricted visible part of the stellar surface might not be rep- 
resentative of the whole star if, for example, it is covered by pre- 
dominantly horizontal magnetic fields and by spots of enhanced 
iron-peak element abundance. Precise spectroscopic and polar- 
ization observations of a Cir over the entire rotation cycle are 
required to investigate the surface magnetic and abundance ge- 
ometry. 
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Fig. 4. Same as in Fig.[T]but for Cr. 
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Fig. 7. Theoretical spectral energy distributions for the model atmosphere parameters T^ff - 7500 K, \ogg -4.1 and different chem- 
ical compositions of aCir. The LLmodels flux spectra, convolved with a FWHM = 10 A Gaussian profile, are presented for our 
final self-consistent model atmosphere with stratification included (black line), chemically-homogeneous solar abundance model 
(red line), and a model computed for BNC08 abundances (blue line). 
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Fig. 10. Comparison of the observed and theoretical Ha and Hy6 line profiles of or Cir Observations are shown by the black line. 
Theoretical calculations with parameters T^s - 7500 K, log g-A-.l are presented for our final self-consistent model atmosphere with 
stratification included {blue line), chemically-homogeneous solar abundance model (green line), and a model computed for lBNCOSi 
abundances (red line). 



